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ABSTRACT 

• Tadpole galaxies, with a bright peripheral clump on a faint tail, are morphological types unusual in 
I I the nearby universe but very common early on. Low mass local tadpoles were identified and studied 

, photometrically in a previous work, which we complete here analyzing their chemical and dynamical 

■ properties. We measure Ha velocity curves of seven local tadpoles, representing 50 % of the initial 

' sample. Five of them show evidence for rotation 70%), and a sixth target hints at it. Often 

, the center of rotation is spatially offset with respect to the tadpole head (3 out of 5 cases). The 

size and velocity dispersion of the heads arc typical of giant HII regions, and three of them yield 
dynamical masses in fair agreement with their stellar masses as inferred from photometry. In four 
OO ', cases the velocity dispersion at the head is reduced with respect to its immediate surroundings. The 

oxygen metallicity estimated from [Nil] A6583/Ha often shows significant spatial variations across 
the galaxies (~0.5 dex), being smallest at the head and larger elsewhere. The resulting chemical 
abundance gradients are opposite to the ones observed in local spirals, but agrees with disk galaxies 
at high redshift. We interpret the metallicity variation as a sign of external gas accretion (cold- 
' flows) onto the head of the tadpole. The galaxies are low metallicity outliers of the mass-metallicity 

relationship. In particular, two of the tadpole heads are extremely metal poor, with a metallicity 
smaller than a tenth of the solar value. These two targets are also very young (ages smaller than 
5 Myr) . All these results combined are consistent with the local tadpole galaxies being disks in early 
Q I stages of assembling, with their star formation sustained by accretion of external metal poor gas. 

Subject headings: galaxies: abundances - galaxies: dwarf - galaxies: evolution - galaxies: formation 
^ I - galaxies: kinematics and dynamics - galaxies: structure 

1. INTRODUCTION cal content. Extremely Metal Poor (XMP) galaxies are 

" > : Tadpole galaxies consist of a large star-forming clump the least evolved objects in the local universe (e.g. 

CN ■ at one end (the head) and a long diffuse region to one Page! et al.l[l992tlKunth fc Ostlinl[200a llzotov fc ThuanI 

^ ] side (the tail). This asymmetric morphology is rather [2QQD- They represent only 1% of the galaxies in a n 

'. common at high redshift but rare in the local universe. arbitrary nearby volume (e.g., |Morales-Luis et aD|201lD, 

• For example, tadpoles constitute 10% of all galaxies but a significant fraction of these chemically primitive ob- 
; larger t han 10 pixels in the Hubble Ultra Deep Field jects turn out to have tadpole or cometary shape (~ 75 %; 

O' (UDF ) (|Elmegreen et 311120071 : IMmegreen fc Elmegreenl iPapaderos et all l2008t iMorales-Luis et all 1201 ID. This 

Cn ■ I2010( ). an d they represent 6% of the U DF galaxies iden- association between low metallicity and tadpole shape 

tified by iStraughn erall (|2006[ ) and [Windhorst et al.l suggests that they are attributes characteristic of very 

(I2006D using aut o mated search algorithms. In contrast, young systems. 

lElmegreen et al.l (|2012bl hereafter Paper I) find only The t adpole structure could h a ve a variety of ori- 

' 0.2% tadpoles amo ng the star-forming lo c al gal axies of gins. [Elmegreen fc Elmegreej (|2010D showed lop- 

j_j ! the Kiso survev bv Mivauchi-Isobe et al.l (|2010( ). This sided ring- like galaxies that would look like tad- 

C3 . decrease suggests the tadpole morphology to represent poles if viewed edge-on. Tadpoles could also result 

a common but transit phase during the assembly of from mergers or close galaxy-galaxy interactions (e.g. 

some galaxies. Since local tadpole galaxies are very low Baldwin ct al. 1982; Corbin et al. 2005; 'Straughn et alj 

mass objects compared to their high redshift analogues |2Q06; Windhorst et al. 2006; D'Onghia et al. 201Qi). 

(jElmegreen et al.ll2012bh . the phase must be alreadv over However, this merger scenario cannot be the universal 

for the local descendants of high redshift tadpoles. pathway to tadpole formation, at least not locally, since 

Another independent observation also suggests that the class of dwarf local tadpole s tend s to be relatively 



o 



the tadpole morphology characterizes a very early phase isolat ed (jCampos-Aguilar et al.l 119931 : iKunth fc Ostlin 

Ca] 



of evolution. From the point of view of their chemi- |200(3| ) and lacks obvious tidal features ([Papadcro s et al 



2008). Another possibility is that the lopsided starburst 

|jos@iac.es, c mt@iac.es, elmegreen@vassar.edu, bge@us.ibm.com, jairoilg mitej from ram com pression by motion through th e 
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stripped gas (e.g.. iChung et al.ll2009D . The scenario m 
which dwarf galaxies are converted in to dwarf spheroidal 
galaxie s in the halos of larger galaxies (|Lin fc Fabeiill983l : 
Ivan den Bcr gh 199 llMa.ver et al.ll2006l[ or in g alaxy clus- 
ters (|BoseUi et airi2008l : Ivan Zee erall 120041 ) could in- 
volve a tadpole phase as the gas and young stars are 
pulled behind. A third possibility is that some tad- 
poles are normal disk galaxies with a large turbulent 
Jeans length for gravitational collapse of their interstellar 
medium. This happens in galaxies that have either small 
rotational velocities or large turbulent motions, and in 
this case only large star f ormi ng clumps can be produ ced 
(|Forster Schreiber et al.ll2006l : lElmegreen et aLlboOQl) . If 
there is only one clump at any particular time, then it 
will appear as a tadpole viewed from the right perspec- 
tive. A fourth possibility is that the head-tail structure 
resul ts from propagating star formation acro ss a disk 
(e.g.. iFranx et alJll997t iPapaderos et al.ll2008| ). Finally, 
as we propose in this paper, the head of the tadpole 
may result from accretion of external flows of pristine 
gas, that penetrate the dark matter halo and hit and 
heat a pre-existing disk. Cosmological simulations pre- 
dict cold-fl ow buildup to be the main mode of galaxy 
formation (|Dekel et al.ll2009l : iGenel et al.ll2012[ ). and the 
incoming gas is expected to form giant cl umps that spi- 
ral in and merge into a central spheroid (jNoguchil 119991 : 
iGenzel et all 120081: lEhnegreen et al.ll2008D . In addition, 
inflow of low-metallicity gas seems to be suggested by 
distorted HI velocities a nd morp hologies of some local 
star-forming galaxies (see lLopez-Sanchcz et al.ii20l3 . and 
references therein). 

It is important to realize that the different formation 
mechanisms mentioned above are not mutually exclusive, 
but exhibit a large degree of overlap between them. For 
example, the cold-flow accretion may be considered as a 
(minor-)merger accretion, and it may also be regarded 
as the interaction of the galaxy with the (fllamcntary) 
intergalactic medium. 

In this context of galaxy formation, we examined a 
sample of fourteen tadpole galaxies in the local uni- 
verse for comparison with high-rcdshift tadpoles (Pa- 
per I). These local tadpoles seem to form a con- 
tinuous sequence with t he U DF tadpoles studied by 
lElmegreen fc ElmegreenI (|2010f ). With regards to their 
photometric properties, local tadpoles occupy the low 
mass end in sequences such as star formation, surface 
density and mass-to-light ratio. In addition, the radial 
intensity profiles of the local tadpoles show an exponen- 
tial decrease at large galactocentric distances, which was 
interpreted as evidence for an underlying disk. The work 
in Paper I is followed up in the present paper. In order 
to determine the dynamical properties and metallicities 
of local tadpoles, we measure Ha spectra along the head- 
tail direction in a representative fraction of the original 
sample. The results reported here show the galaxies to be 
rotating structures with an unexpected metallicity pat- 
tern, which may shed light into the nature of the tadpole 
morphology. 

The paper is organized as follows: Sect. [2] describes 
the observations, the main steps of the reduction, and 
the main properties of the resulting spectra. Section [3] 
details the approximations used along the work to derive 
physical parameters from spectra. Section |4] analyzes 
the dynamical properties of the galaxies, their rotation 



curves (RCs) and linewidths. These are used to infer dy- 
namical masses, which are then compared with the stellar 
masses derived in Paper I (Sect.[6l). Section [5] studies the 
light profiles both in continuum and Ha emission. The 
variation of the oxygen abundance along the galaxy is 
measured in Sect. [T] Notes on individual galaxies are 
given in Sect. [5] The results of our observation in the 
context of tadpoles as dynamically young disks are dis- 
cussed in Sect. |9l 

2. OBSERVATIONS AND DATA REDUCTION 

The sample of local tadpole galaxies selected in Pa- 
per I comes from the Kiso s urvey of UV intense galaxies 
(|Mivauchi-Isobe et al.ll20Tol) . Among those galaxies la- 
belled in the Kiso catalog as clumpy or having conspicu- 
ous HII regi ons, we visually inspect ed 158 with images in 
SDSS-DR7 (jAbazaiian et al.ll2009[ ). Then tadpoles were 
subjectively selected as those objects with lopsided light 
distribution, so that the brightest clump is far to one end 
and the rest of the galaxy is mostly featureless. This in- 
spection rendered only 13 targets, which we completed 
with one additional galaxy with similar features from the 
University of M ichigan survey of emission-line objects 
(jMacAlpine fc L ewis 197^. Half of the sample in Pa- 
per I was randomly selected for the follow-up work pre- 
sented in this paper. The observing program was aimed 
at determining the dynamical properties of the galaxies 
based on long slit spectroscopic observations along the 
well-defined head-tail direction (see Fig. [Ij. The project 
was planned for the spectrograph IDS of the 2.5 m Isaac 
Newton Telescope (I NT) at the Observat orio del Roque 
de Los Muchachos (|Laing fc JonesI I1985D . It provides 
adequate spectral resolution to determine centroids of 
emission lines with an accuracy better than lOkms^^, 
which suffices to characterize RCs even for dwarf galax- 
ies. The original project could not be completed due to 
poor weather conditions, so it was finished using service 
time of the 2.5m Nordic Optic al Telescope (NO T), op- 
erated in the same observatory ()Andersenlll985l ). 

The observing logbook is summarized in Table [TJ 
The used IDS@INT setup includes a grating with 
1800 lines mm^^, which provides a scale in the spec- 
trograph focal plane of 0.32Apix~^ equivalent to 
0'.'44pix^^. The spectral resolution is set by the width of 
the 1" slit, which corresponds to 2.2 pix, or 0.73 A, or 33 
kms^^ in Ha. The camera covers some 700 A around the 
central wavelength, which was tuned to Ha, and which 
automatically included the line [NII]A6583 used in our 
metallicity estimates (Sect. The tadpoles KIS03193, 
KIS03867, KIS05149 and KIS08466 were observed with 
this setup (Table [T|). Figure [5] shows one of them as an 
example. The image contains the spectral region around 
Ha. The reduction procedure included the standard bias 
and flatfield corrections, cosmic ray elimination, as well 
as removal of sky emission lines. The typical signal-to- 
noise (S/N) ratio in the peak Ha emission exceeds 500, 
and it decreases down to one in the outskirts of the galaxy 
(Fig. O bottom panel). This S/N is achieved after inte- 
grating for some one and a half hours (Table [1]). Seeing 
was of the order of 1" (Table [T]). 

The other three galaxies, KIS05639, KIS06669, and 
KIS06877 were observed with the NOT telescope with an 
observational setup equivalent to the previous one. We 
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Fig. 1. — Images of the observed galaxies with the red line showing the orientation of the slit and with the arrow indicating the sense 
of increasing distance along the slit used in the paper. The images have been taken from the SDSS database, and are displayed with an 
inverted color palette so that the background sky looks white, and the intrinsically blue galaxies appear reddish. The horizontal scales on 
the upper left corner of the panels correspond to 5"or 10"as indicated. 



used the ALFOSC spectrograph which, together with a 
high resolution grism, provided a scale of 0.26Apix~^x 
0'.'19pix~^ at Ha. This pixel size over-samples the spec- 
trum, whose resolution was set by the CK'Q wide slit. In 
order to match the pixel size to the actual resolution, 
the original data were re-binned 2x2, yielding the fi- 
nal resolution indicated in Table [TJ The data reduction 
followed the standard process mentioned above, with a 
small difference worthwhile mentioning. The standard 
procedure includes using images of the sky taken during 
twilight to characterize (and then to correct for) large- 
scale illumination gradients in the images. Such sky-flat 
images corresponding to KIS06877 were not available. 
For reduction we used those from the previous night, 
but this rendered spectra with small residual sky lines. 



To estimate the uncertainties introduced by this prob- 
lem, the reduction was independently repeated by two of 
us (JSA and JMA). Thus different spectral regions and 
interpolations were used to construct the flatfield and to 
subtract the sky background. The two reductions are in 
good agreement and provide consistent results. Table [2] 
contains the final physical parameters inferred from the 
two parallel reductions, and they agree to the point that 
their small differences do not affect the results in the 
work. Seeing during observations was sub-arcsec. The 
NOT spectra have S/N similar to the INT spectra; the 
spectral resolution is slightly worse, and the angular res- 
olution slightly better. 

Determining the actual spectral resolution is central 
for the proper interpretation of the observed linewidths. 



4 



6700 - 
6690 ' 
6680 ^ 
6670 ^ 
6660 
6650 ' 
6640 ^_ 




-10 



10 



20 



-20 

Tarcsecl 



20 



40 



kiso8466 



600 — 



400 — 




200 — 



6640 



6660 6680 
Wavelength [A] 



6700 



6720 



Fig. 2. — Example of a reduced spectrum from IDS@INT. The image shows a negative of the spectrum after reduction - x-axis represents 
position and y-axis wavelength. Positions are given in arcsec from the galaxy head, whereas wavelengths arc in A. The wavelength integrated 
flux (spectral flux) is shown in the top right corner. Spectra corresponding to various positions are shown in the bottom panel, with the 
positions indicated as vertical dotted lines in the two upper plots. The example corresponds to KIS05149, but it is representative of the 
full set. The largest line is Ha whereas the second largest is [NIIJA6583, used in our metallicity estimates. 



As usual, it was determined from the width of the telluric 
lines in the spectra since their intrinsic widths are just 
a fe w kms~^ and thus negligible (e.g.. iBalthasar et alj 
Il982i) . Gaussian functions were fitted to the lines, 
and their mean Full Width Half Maximum (FWHM) 
are quoted in Table [1] as the spectral resolution, with 
the error bars representing the rms fluctuations along 
the field of view. We average over all the targets ob- 
served with the same instrumental configuration to es- 
timate the resolution corresponding to that particular 
configuration. The weighted average for the IDS@INT 
spectra turns out to be 25.4 ± 0.6 kms~^, whereas in 
the case of ALFOSC@NOT the measured resolution is 
54.4 ± 3.9 kms~^. The error bars represent the weighted 
rms fluctuations. 

3. DETERMINATION OF PHYSICAL 
PARAMETERS 

Velocities, masses, abundances and other physical pa- 
rameters are determined from the spectra. This section 
explains how these parameters are computed, including 
the underlying hypotheses and their uncertainties. 

Bulk velocities arc measured from the displacement of 
Ha. Wc compute the displacement both as the bary cen- 



ter of the emission line, and as the center of a Gaussian 
function fltted to the profile. Errors are estimated from 
the S/N measured in the continuum and then propa- 
gated to the centroids (e.g.,[N'Iartin 1971i Sect. 5.3). The 
FWHM of the profiles are also measured directly from 
the profile and from the Gaussian fit. Their errors are 
also inferred from the noise measured in the continuum 
by error propagation. 

The variation of the velocity with position along the 
slit is usually referred to as the velocity curve. The ve- 
locity curves of our targets can be ascribed to rotation. 
In order to characterize such rotation, we fit an analytic 
RC to the velocity curve. We wanted a function that is 
simple but produces a good m atch to the observa tions. 
The universal RC advocated bv lSalucci et al.l (j2007| ) does 
a good job representing the observed variation. Specifi- 
cally, we adopted the dark matter component of the uni- 
versal RC given by, 

u{d) = uo + u, J-'^; (1) 

where U{d) is the velocity observed at a distance d, and 
Uo,Ui, A, and a-re free parameters to be determined 
from a non- linear fit. The curve is simple - do yields a 
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TABLE 1 

Observing Logbook 



Name^ 


RA DEC 


t-cxp 


Date 


Instrument 


Slit 


Seeing'^ 


Pixel 


Scale 


Tclluric<^ 






[s] 


2012 










[pc/"] 


[kins^^] 


KIS03193 


08:56:08 +39:52:09 


6000 


Feb 10 


IDSQINT 


I'.'O 


I'.'O 


0.32AxO'.'44 


124 


26.0+1.9 


KIS03867 


09:40:13 +29:35:30 


6000 


Feb 10 


IDSQINT 


I'.'O 


I'.'O 


0.321x0^44 


36 


25.0+1.0 


KIS05149 


11:16:08 +23:29:16 


6000 


Feb 10 


IDSQINT 


I'.'O 


I'.'O 


0.32Ax0'.'44 


834 


28.2+1.5 


KIS05639 


11:41:07 +32:25:37 


5000 


May 29 


ALFOSC@NOT 


o'.'g 


0'.'8 


0.511x0^38 


119 


53.9+5.7 


KIS06669 


12:31:50 +27:23:13 


5000 


May 29 


ALFOSC@NOT 


O'.'Q 


0'.'8 


0.5lAx0y38 


299 


55.6+8.7 


KIS06877 


12:46:11 +26:15:01 


5000 


May 31 


ALFOSC@NOT 


O'.'Q 


0'.'7 


0.5lAx0('38 


128 


54.3+6.6 


KIS08466 


16:03:27 +19:09:46 


4800 


Feb 10 


IDSOINT 


I'.'O 


I'.'O 


0.32Ax0^'44 


313 


24.9+0.7 



^ Named as in lElmegreen et al.l I l2012bl) . Paper I. 
^ Total exposure time. 

^ Mean value for the night from RoboDIMM@WHT; see 'http : / /cat server . ing. iac . es/robodinun/j 
^ FWHM of the observed telluric lines, used as proxy for spectral resolution. 



center for the RC, Ui gives its amplitude (t/ — > Uq±Ui 
when d — >■ ±(X)), and A provides the spatial scale for the 
central gradient. We stress that the hypotheses behind 
the analytic expression ([T]) are of little importance in our 
context. The formula is used here because it provides a 
good representation of the observations, and so provides 
a smoothed version of the velocity curve. Moreover, it 
allows us to determine the center of rotation do- The 
center thus determined does not depend so much on the 
actual expression used to parameterize the RC, but on 
the fact that the curvature changes sign at the center of 
rotation, i.e., 

f/"(rfo) - 0, (2) 

with U" the second derivative of U . This is a general 
property that any anti-symmetric RC satisfies. 

Dynamical masses are estimated from RCs and 
lincwidths. Assuming the mass distribution to be spheri- 
cally symmetric (e.g., a bulge or a dark matter halo), the 
circular velocity Uc that balances the gravitational pull 
depends only on the mass enclosed within the radius p, 
M{p), 

Ul{p) = GM{p)/p, (3) 

with G the gravitational constant. We assume that the 
measured macroscopic velocities are the circular veloci- 
ties affected by the inclination i of the galaxy plane, i.e., 
U = Uc sini. We also assume the spectrograph slit to be 
oriented along the galaxy major axis, implying p = d—do. 
Then the inner mass of the galaxy up to distance d has 
the usual expression, 

M (d) sin^ i^{d- da) U^{d)/G, (4) 

which using astrophysical units turns out to be 

M{d) sin^ i = (2.33 • IO^Mq) {d - do) U{d)^, (5) 

with d measured in kpc and U in kms~^. Equa- 
tion ([5|) is an approximation; however, including more 
realism in the m ass distribution (e.g., pressure support, 
Dalcanton fc Stilp. ,20 1 0i: or non-spherical components, 
Salucci &: Persid Il997| ) would only modify the scaling 
factor as a correction of order one. Masses of individ- 
ual clumps are inferred from linewidths assuming virial 
equilibrium. If the isotropic velocity distribution that 
balances the clump gravity has a dispersion a, then 



(|Bosch et al.|[2009l . Sect. 4.2) 

M = (1.20 • IO^Mq) Re FWHM^ (6) 

where 

FWHM = 2\/21n2cr. (7) 

The symbol stands for the half-light radius of the 
clump, so that Eq. ([6|) is equivalent to Eq. ([5]) after in- 
cluding the appropriate scaling factors. As in the case 
of Eq. ([5]) , distances are measured in kpc and linewidths 
in kms^^. The linewidth in Eq. ([5]) is not the observed 
width FWHMo, but the width corrected for instrumen- 
tal spectral resolution FWHMi, thermal motions in the 
nebula FWHMt, and natural width of Ha FWHMn, 

FWHM^ FWHM^ - FWHM? - FWHM^ - FWHM^, 

, , (8)_ 

see, e.g., iTerlevich fc Melnickl (|1981| ): iMelnick et all 
([19991) . We use the measured widths of the telluric lines 
as proxy for instrumental broadening (Table [T]). The 
thermal broadening is assumed be the same for all galax- 
ies at all positions. Its value has been set to a representa- 
tive round number FWHMt = 25kms~"'^, which approx- 
imately corresponds to H atoms at 14000 K, a temper- 
ature typical of HII regions. The actual FWHMt is of 
secondary importance since the range of possible values 
is significantly smaller than the FWHM resulting from 
Eq. (|8]E|. Finally, the natural width of Ha is of the or- 
der of 7kms~^ (e.g-. lRozas et al.ll2006[) . The radius Re 
used to estimate dynamical masses from line widths using 
Eq. © is computed fitting a 1-D Gaussian to the light 
distribution across the galaxy head. The observed light 
profile is assumed to represent a 1-D cut across the center 
of a 2D-Gaussian, which readily provides the half-light 
radius from the width of the fitted Gaussian - in a 2D 
Gaussian, the FWHM is twice the half-light radius. The 
observed radii are corrected for seeing using a formula 

^ FWHMt varies only from 19 to and 26 kms~^ for temperatures 
between 8000 and 15000 K. Moreover, we did the exercise of calcu- 
lating masses also with the extreme value of FWHMt=0, to check 
that this assumption does not modify the conclusions in the work. 
Setting FWHMt =0 is equivalent to including thermal motions as 
part of the virial equilibrium represented by Eq. The exercise 
implies that we do not have to worry about whether the thermal 
motions contribute or not to the virial equilibrium. Including or not 
thermal motions as part of the kinetic energy does not significantly 
modify the masses estimated in the work. 
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similar to Eq. (jS]), namely, 



eO 



{S/2f 



(9) 



where i?eo stands the measured effective radius and S 
represents the seeing, i.e., the FWHM of the seeing disk 
as given in Table [TJ 

The metallicity of the gas is commonly inferred by com- 
bining emission-line fluxes of several atomic species to 
derive their relative abu ndances (e.g.. lPagel fc Edmundsl 
119811: IOsterbroclj|1989D . This approach is the so-called 
direct method or temperature-based method, and it is 
to be preferred whenever possible. However, it involves 
measuring fluxes of lines spread throughout the visiblc- 
IR spectrum, so is expensive observationally. Fortu- 
natel y, we have alternatives called strong-li ne methods 
fe.g.. lShi et aII[200llKe^lev fc Ellisonl[200l . where the 
metallicity is estimated empirically by relating the ra- 
tio of a few selected line fluxes with the abundance of a 
particularly relevant metal ( typica lly oxygen). The one 
proposed bv lDenicolo et all (|2002D turns out to be ideal 
in our case, when only the spectral region around Ha is 
available. It yields the oxygen abundance from the ratio 
of [NIIJA6583 to Ha, and [NIIJA6583 automatically ap- 
pears in the s pectra next to Ha (see Fig: . p). \y e use the 
calibration bv iPerez-Montero fc Continil ()2009t) . 

12-Hlog(0/H) = 9.07 0.79 log([NII]A6583/Ha), (10) 

particularly suited f or low metallicity targets (c.f., 
IPettini k, Pagell I2004D . Equation ([TO]) provides the O 
metallicity from the flux in a N line, therefore, it may be 
biased in objects having unusual N/0. In order to discard 
this potential bias in our O metallicities, we also estimate 
the ratio of N to O using the sulfur lines [SIIJA6717 and 
[SIIJA6731 present in most of our spectra. We use the 
calibration 

log(N/0) = 1.261og([NII]A6583/[SII]A6717,6731)-0.86, 

, , (11) 

also by iPerez-Montero fc Continil (I2009D . 

The spectral flux of the galaxy is computed by inte- 
grating the observed spectra between their two extreme 
wavelengths (from 6242 A to 6935 A for INT spectra, and 
from 6368 A to 6840 A for NOT spectra). The flux in 
Ha results from integration of the emission line profile 
around its maximum, once the underlying continuum was 
removed. The integration includes a 10 A wide region 
around the peak emission, whereas the continuum was 
obtained by fitting a linear function to two continuum 
windows, 10 A wide, outside the line. The Ha equiva- 
lent width (EW) is inferred from the Ha flux dividing by 
this continuum. Continuum and emission lines are com- 
bined in the spectral flux for simplicity, and it suffices to 
compare the limited extension of the emission line region 
inferred from Ha with the rest of the galaxy. 

4. VELOCITY CURVES AND LINEWIDTHS 

Figure 13] shows the velocity curves of the tadpole galax- 
ies included in our study. Five out of the seven targets 
show velocity gradients interpreted as rotation. The fig- 
ure includes the best fit to the analytic RC in Eq. ([T]), 
which does a good job reproducing the observations. We 
use for fitting the portion of the velocity curve indicated 
in red in Fig. |3l which includes all positions but the 



extremes with large error bars or obvious distortions. 
Sometimes the interpretation of the velocity curve as 
a RC is obvious (e.g., KIS08466), but other times the 
curve looks more like a perturbed RC (e.g., KIS05639). 
KIS03193 and KIS03867 have a rather flat velocity curve, 
and therefore no obvious rotation. However, one of them, 
KIS03867, shows a systematic line shift of the order 
10-20 kms~^ between the two extremes of the galaxy 
(Fig. [3]). The amplitude is of the order of the error bars, 
but the displacement is in the raw data as judged by in- 
spection of the individual Ha profiles. Table [2] contains 
the amplitudes of the RCs as assigned by the fit (i.e., Ui 
in Eq. [T]). Figure [3] represents velocities obtained from 
the baryccnter of Ha (Sect. |3]). The velocities from the 
Gaussian fit arc not shown because the two estimates dif- 
fer only by a few kms~^, a difference always smaller than 
the error bar assigned to each velocity measurement. 

Table also contains the center of rotation do , which 
often differs from zero, i.e., from the position of the tad- 
pole head. As argued in Sect. [3l its estimate is fairly 
robust since it comes from the point where the curva- 
ture of the RC changes sign - accordingly, it has the 
small formal error bars provided in Table [51 Figure [3] 
shows that three out of the five rotating galaxies, ex- 
phcitly KIS06669, KIS06877 and KIS08466, have their 
center of rotation displaced with respect to the head by 
more than the 1" uncertainty introduced by seeing (see 
Table [1}. Considering that the rotation center points out 
the center of the galaxy, our results indicate that often 
the star-forming region at the head is displaced with re- 
spect to the center of the galaxy. This suggests that the 
head is not a bulge-like central spheroid. The fact that 
the center of rotation is sometimes displaced from the 
heads may be easier to appreciate in Fig. 21 

Figure |S1 shows the FWHM of Ha as inferred directly 
from the emission line profile (the solid lines) and from 
the Gaussian fit (the dashed lines). A significant part of 
the observed linewidths is due to instrumental broaden- 
ing. The widths corrected for instrumental effects and 
thermal motions are given as dotted lines in the figure. 
We use the directly inferred widths to estimate the in- 
trinsic widths, but using them or the widths from the 
Gaussian fit render similar results since the two mea- 
surements are in very good agreement (Fig. [5]). First 
note that the widths of all galaxies, large and small, are 
of the order of 20-70 kms"^. These widths are typi- 
cal of giant HII regions such as 30 Doradus rather than 
the widths of the HII regions observed in large spirals 
(e.g.. lO'deU fc Townslevl[l988tlMufioz-fiIfiOT]ll99l . Fig- 
ure ini shows the variation of the tadpole head size as 
a function of its velocity dispersion. The two quanti- 
ties are known to be correlated in giant HII regions, so 
that the larger the dispersion the larger the size (e.g., 
lMuiioz-Tufioa[T99lll\ientes-Masip et al.ll2000( ). We find 
the tadpole heads to f ollow s u ch trend as charac terized by 
iTerlevich fc Me lnick" (TgSllV IRov et al.l (fT986l) . or more 
recentlv bv IWisnioski et al.l (|2012D (see Fig. [Q. In the 
case of KIS06877, the error bars of the line width at the 
head are so large that the width is actually an upper 
limit, but even with this caveat in mind, the measure- 
ment is consistent with its head being a giant HII region 
(Fig.®. 

The FWHM also varies along the galaxies, and the 
fluctuations are correlated neither with the spectral flux 
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TABLE 2 
Physical parameters 



Name^ 


Distance^ 


FWHM'' 


logMh D'= 




logMh Ph<= 




do*5 


log M D'^ 


log M Ph' 


12+log(0/H)j 




[Mpc] 


[kms^-*-] 




[kpc] 




[kms 


[kpc] 


[Me] 


[Mq] 




KIS03193 


25.5 


39.3 ±2.0 


7.61 ± 0.04 


0.221 ± 0.009 


6.4 ±0.1 








7.3 ±0.4 


7.89 ±0.31 


KIS03867 


7.4 


35.4 ± 2.2 


7.22 ±0.04 


0.109 ± 0.003 


5.2 ±0.3 








6.8 ±0.4 


8.03 ±0.26 


KIS05149 


172 


70.1 ±3.0 


9.03 ±0.03 


1.825 ± 0.064 


8.9 ±0.1 


168 ± 92 


0.3 ±3.0 


10.8 ±0.6 


9.6 ±0.4 


8.42 ±0.06 


KIS05639 


24.5 


46 ± 11 


7.22 ±0.15 


0.065 ± 0.004 


6.7 ± 0.2 


34.7 ± 6.2 


0.180 ±0.064 


8.2 ±0.2 


7.7 ± 0.4 


7.48 ±0.04 


KIS06669 


61.6 


47 ± 11 


7.73 ±0.14 


0.206 ± 0.011 


7.0 ± 0.1 


71.4 ± 4.6 


1.24 ± 0.16 


9.8 ±0.1 


8.9 ±0.4 


8.11 ± 0.12 


KIS06877 


26.3 


22 ± 23 


6.19 ±0.66 


0.028 ± 0.004 


6.3 ±0.1 


19.0 ± 1.7 


0.246 ±0.020 


7.7 ± 0.1 


7.3 ±0.4 


7.40 ± 0.09 


Kiso6877'^ 


26.3 


21 ± 24 


6.30 ±0.70 


0.038 ± 0.003 


6.3 ±0.1 


19.0 ± 1.6 


0.242 ±0.020 


7.6 ±0.1 


7.3 ±0.4 


7.41 ±0.07 


KIS08466 


64.5 


61.0 ±2.9 


8.19 ±0.04 


0.345 ± 0.021 


8.1 ± 0.1 


117 ± 16 


2.95 ±0.36 


10.0 ±0.1 


9.2 ±0.3 


8.14 ± 0.05 


" FromlElmeereen et alj (I2012bl), Paper 


I. 

















Ha width at the tadpole head. 
Dynamical mass of the head — from Eq. JHl . 
^ Half-light radius of the head. 

^ Photometric mass of the head, from Paper I, Table 3. 

^ Amplitude of the rotation curve — see Eq. fll . 

^ Center of rotation relative to the tadpole head - sec Eq. 0. 

^ Dynamical mass of the galaxy - sec text for its computation. 

^ Photometric mass of the galaxy, from Paper I, Table 3. 

J Metallicity at the head - from Eq. ilOt . 

^ Second independent reduction. 

nor with the Ha flux (c.f. Fig. |4]and[5]). If anything, 
there is a tendency for the tadpole heads to coincide with 
local minima of lincwidth (see KIS05149, KIS05639, and 
KIS06877 in Fig. [5]). Similarly, the center of rotation 
of the rotating galaxies docs not seem to be associated 
with extremes of the FWHM curve. Interpreting 
the decrease of linewidth associated with the head is not 
straightforward. One may try to explain its origin in 
the context of forming stars in a highly turbulent galaxy. 
Only where the turbulence is low enough the conditions 
that trigger star formation are met, and this preference 
for low turbulence regions is what we detect. 

The square of the ratio between velocity dispersion and 
rotati on gives the Jeans length relat ive to the galaxy size 
(e.g., iBournaud fc Elmegreen|[2"009[ ). In other words, it 
provides the relative size of the clumps to be produced 
by gravitational instability. Using Eq. ([7]) and the data 
in Table [21 the ratio dispersion to rotation a/Ui turns 
out to be between 0.2 and 0.6 at the tadpole heads, 
which yields expected clump sizes between 0.03% and 
30 % galaxy radii. The relative size increases with de- 
creasing galaxy mass, and it may be a coincidence, but 
the largest ratios correspond to the two XMP galaxies 
in the sample (see Sect. [7]). In some cases the predicted 
clumps have sizes comparable to those of the observed 
heads. The ratio a/Ui also measures whether the object 
is supported by random motions or rotation. Our disks 
have ratios similar to the turbulence-su pported clumpy 
galax ies observed at high redshift (e.g., iCeverino et al.l 

5. FLUXES AND EQUIVALENT WIDTHS 

Figure 13] shows the variation across the galaxies of the 
spectral flux and the Ha flux (the solid lines and the 
dashed lines, respectively). The position where the two 
fluxes are largest coincides. (As we already pointed out, 
the origin of distances used along the paper has been 
set at the maximum of the spectral flux distribution.) 
Note also the obvious lopsidedness of all light curves, as 
expected in tadpole galaxies. Another important fea- 
ture is the extension of the Ha emission as compared to 



the spectral flux emission which includes Ha plus contin- 
uum. Ha is more concentrated, implying that the region 
producing the total emission extends further away from 
the star-forming regions. In other words, it shows the 
existence of an underlying galaxy, with old stellar pop- 
ulations that contribute to the spectral flux but not to 
Ha. 

The dotted line in Fig. H] represents a Gaussian plus 
a constant term fitted to the Ha flux around the head. 
The fits are good, and they allow us to assign a radius 
to the head. Table [2] lists the effective radii of the heads 
Re, defined as the radius enclosing half of their light (see 
Sect. [3]). They span a wide range of values from 1.8 kpc to 
50 pc, refiecting the wide range of intrinsic galaxy sizes. 

Figure [7] shows the variation of Ha EWs across the 
galaxies. The largest EWs tend to coincide with the max- 
imum fiux (Fig. |4|), although not always (kiso3193 in 
Fig. [7|). In general, EWs are fairly moderate, with max- 
ima <300A. The two exceptions correspond to the two 
extremely metal poor targets to be described in Sect. [7]- 
KIS05639 and KIS06877. Their large Ha EW implies the 
extreme youth of the star-forming regions at the galaxy 
head. We know from modeling that the Ha EW of an HII 
region must be smaller than some 3000 A, and it drops 
down very quickly so that a coeval starburst reaches 
EW~200A in just lOMyr (e.g.. iLeitherer et al.l[T99l 
Fig. 83). There is also a dependence of the EW on metal- 
licity, but age is by far the dominant factor. The large 
observed Ha EWs correspond to ages of only a few Myr, 
and these ages arc upper limits to the stars responsible 
for the ionization since the (old) underlying galaxy pro- 
duces continuum emission that reduces the observed EW. 
If KIS05639 and KIS06877 are as young as we infer from 
their Ha EWs, one expects to find Wolf-Rayet (WR) 
star features in the spectra, which are characteristics 
of ex tremely young starbursts (< 5 Myr; e.g.. iCrowtherl 
|2007| ). These features are distinctive broad bumps at 
4600-4680 A and 5650-5 800 A fe.g.. lSchaerer et al.l[l99l 
iBrinchmann et afl l2008l) . Unfortunately, the WR fea- 
tures lie outside our observed spectral range. However, 
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Fig. 3. — Velocity curves of the seven tadpole galaxies. The abscissae represent distances along the major axes of the galaxies from 
the position of the tadpole head; i.e., the brightest point on the galaxy. The range of distances differs for the different targets, but the 
horizontal bar in each panel gives a common length scale corresponding to 1 kpc. The points with error bars show the observations whereas 
the thick solid line represents the best fit of the observed points to the analytic RC. The part of the RC shown in red indicates the portion of 
the velocity curve used for fitting. The thin horizontal and vertical lines indicate the systemic velocity and the center of rotation obtained 
from the fit, respectively. The little arrows on top of each panel also indicates the center of rotation. The zero of the velocity scale is set 
by velocity of the spatially integrated spectrum of the galaxy. Positive velocities are redshifts. The ordering of the galaxies is identical to 
that in Fig. [T] 



our ga laxies also have SPSS spectra, which coyer a wider 
range (|Stoughton et al.ll2002tlAbazaiian et"ani2009D . We 
inspected them for WR features but we did not see any 
(see Fig.H]). Moreover, our targets were not found in the 
sys tematic search for WR g alaxies in SDSS carried out 
bv iBrinchmann et all (|2008l) . KIS06877 does not show 
WR features, most probably because the SDSS spec- 
trum was taken away from the tadpole head, in a region 
which is not particularly young. KIS05639 does not show 
the broad WR features either. Instead, its spectrum 
contains narrow high excitation emission lines including 
HeIIA4686 (Fig. [8l). These lines are supposed to be ex- 
cited only by the hard UV-radiation of WR stars, which 
makes interpreting spectr a with HeII A4686 but without 
WR bumps puzzling (see. iShirazi fc Brinchmann. 20121 ). 
However, there is a significant number of star-forming 
galaxies without WR features showing HeIIA4686 emis- 
sion. The reason is unknown, but these galaxie s are usu- 
ally metal-poor (jShirazi fc BrinchmannI [2012|) . so that 
the presence of high excitation narrow lines seems to re- 
flect the extreme youth of a metal-poor starburst. Ac- 



cording to IShirazi fc BrinchmannI ()2012[ ). the stellar pop- 
ulations at very low metallicities can have much higher 
temperatures than is currently expected in models. Then 
even main sequence O stars may excite HeIIA4686. Al- 
ternatively, in low metallicity environments the winds of 
WR stars are weak, and so, optically thin in the He-j- 
continuum, allowing the ionizing radiation to escape cre- 
ating an Helll region respon sible for the observed emis- 
sion (see lKehrig et al.]l201lL and references therein). In 
short, even though KIS05639 lacks WR bumps, the pres- 
ence of high excitation lines such as HeIIA4686 reinforces 
our conclusion that its head contains an extremely young 
metal-poor starburst. 

6. DYNAMICAL MASSES 

Figure [9] shows the relationship between the dynami- 
cal mass and the photometric mass for the tadpole heads 
and the full galaxies. The photometric masses are from 
Paper I, whereas the dynamical masses come from apply- 
ing Eqs. (O and (jB]). In particular, Eq. ([5]) has been inte- 
grated until the last point used to compute the RC. The 
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Fig. 4. — Spectral flux (the solid lines) and Ho flux (the dashed lines). Note the obvious lopsidedness of the light distributions, as 
expected in tadpole galaxies. The origin of distances has been set as the maximum of the spectral flux distribution. The dotted line 
represents a Gaussian fitted to the Ha flux around the head. The horizontal bar in each panel gives a common length scale corresponding 
to 1 kpc. Formal error bars for photometry are not included since they are negligible small. The thin vertical solid lines indicate the center 
of rotation obtained from the RC fit shown in Fig.fS] The ordering of the galaxies and the range of abscissae are identical to those in Fig. [3] 



actual values used for plotting are listed in Table [21 Note 
that all dynamical masses are larger than the photomet- 
ric masses. The difference is less important in the heads 
(except for the case of kiso3867, discussed in Sect. [S]). 
Three of the heads have almost identical dynamical and 
photometric masses, so in these cases there appears to be 
no dark matter in the heads (kiso5149, kiso6877, and 
KIS08466). Keep in mind that the dynamical masses of 
the galaxies derived from RCs are actually lower limits. 
First, there is a sin^ z factor in Eq. ([5]). It is probably 
unimportant since our galaxies are elongated suggesting 
large inclinations and so sin^ i ~ 1. Second, and more 
critical, is the fact that the RC has been integrated only 
to the largest radii having velocities. Since the dynamical 
mass of the galaxies thus derived exceeds the photomet- 
ric mass, we can conclude that galaxies are objects with 
significant amounts of non-stellar matter. 

If the star-forming regions at the head of the tadpole 
were self-gravitating, one would expect them to hold 
some degree of internal rotation. If this rotation sig- 
nificantly differs from the galaxy rotation pattern, and 
if the head is massive enough, then the rotation of the 
head could perturb the rotation curve of the galaxy pro- 



ducing a noticeable distortion (e.g.. Ilmmeli et aT]|2004[) . 
We examined the observed RCs for such signals and did 
not find them, except perhaps in the case of KIS06877, 
discussed in Sect. [S] 

7. METALLICITIES 

The mctallicity was estimated as explained in Sect. |3l 
using the ratio [NIIJA6583 to Ha. Figure [TUl shows the 
variation across the galaxies of the oxygen abundance, in- 
cluding their error bars. (Points with errors larger than 
1 dex have being excluded.) The galaxies tend to have 
sub-solar metallicity (the thick hori zontal line marks the 
solar oxygen abundance given by lAsplund et all [2009I . 
124-log(0/H)Q = 8.69 ±0.05). The galaxies also present 
significant abundance gradients, with the lowest abun- 
dances tending to coincide with the largest Ha emis- 
sions (e.g., KIS06669 and KIS06877 in Fig. [TOl keep- 
ing in mind that the vertical dotted lines mark the po- 
sition of the peak Ha fluxes). We also note that two 
targets, KIS05639 and KIS06877, have metallicities well 
below one-tenth the solar value, therefore, they belong to 
the select ed club of XMP galaxies fe.g.. lKunth fc OstlinI 
[2000; .Guseva et al.|[200l . They are really rare objects: 
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Fig. 5. — FWHM of Ha inferred directly from the emission line profile (the solid lines), from the Gaussian fit (the dashed lines), and 
after correction for instrumental, thermal and natural broadenings (the dotted lines). Formal error bars for the direct measurement are 
represented too. The ordering of the panels and the abscissae are identical to those in Fig. [3] 

one out of a thousand galaxies in th e local universe ac- 
cording to iMorales-Luis et al.l (1201 If ). Therefore the fact 
that we observe two in a sample of seven cannot be 
a coincidence. It is known that a significant fraction 
of XMP galaxies turn out to be cometar y or tadpole 
(|Papaderos et al.l 120081 iMorales-Luis et al.l (20H"l ). Here 
we find that the reverse holds too, i.e., that tadpole 
galaxies have a significant probability of being XMP. As 
we discuss in Sect. |9l this fact supports the idea that the 
tadpole morphology is a sign of dynamical youth, as the 
low metallicity is a sign of being chemically young. 

The observed gradients in metallicity are one of the 
central results of this work and, therefore, deserve a 
separate discussion. Our abundance determinations are 
based on N2=[NII]A6583/Ha rather than on the direct 
meth od, and this ma y be a source of systematic error 
(e.g., iShi et all l2005f ). The spatial gradients in metal- 
licity may be artificially due to gradients in excitation. 
As Morales-Luis et al. (2013, in preparation) discuss, the 
excitation and (to a lesser extent) N/0 change N2 at 
12 -I- log(0/H) ~ 7.5. The higher the excitation the 
smaller N2, and the excitation is expected to change 
with time since the number of ionizing photons drops 
down quickly in young starbursts (e.g., iLeitherer et al.l 
119991 Fig. 77). However, the bias the excitation pro- 
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Fig. 6. — Half-light radius of the tadpole head versus Ha 
linewidth. The velocity dispersions are usually parameterized in 
terms of cr = FWHM/(2V2 In 2), and its scale is given in the 
abscissae on top of the figure. The observed points follow the 
trends observed in giant HII regions - solid, dashed, and the dotted 
lines correspond to the em p irical calibration of th e relationship by 
ITerleyich fc Melnicld lITOSlli IKov et al.l | [T98g) . and lWisnioski et al.l 
II2012I) . respectively. 
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Fig. 7. — Ha EW in A. EWs obtained directly from the line profile (the solid lines) and through a Gaussian fit (the dashed lines) are 
included, although they arc difficult to distinguish. The range of distances differs for the different targets, but the horizontal bar in each 
panel gives a common length scale corresponding to 1 kpc. Only points whore the EW is larger than its formal error bar are included. The 
ordering of the galaxies and the abscissae are identical to those in Fig. [3] 
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Fig. 8. — SDSS spectra of our two youngest targets in the wave- 
length region around one of the WR bumps (4600-4680 A). The 
bump is missing, but one of the targets shows a number of high- 
excitation lines including HeIIA4686. The solid slanted line repre- 
sents the continuum and is included to reinforced the lack of any 
bump above noise. 



Fig. 9. — Dynamical mass versus photometric mass for the heads 
of the galaxies (red square symbols) and the full galaxies (blue as- 
terisks). The slanted solid line shows where the two masses are 
equal. The observed dynamical masses usually exceed the pho- 
tometric masses. The heads of three tadpoles he on the equality 
lino within the errors, meaning that their dynamical masses are 
consistent with their photometric masses. 



duces is much too small to acc ount for the > 0.5 dex 
gradients we detect (Fig. [10]). iSanchez Almeida et alj 



(|2009l Appendix A) studied the difference between the 
oxygen abundance derived from the direct method and 



12 



from N2 in a large set starburst galaxies with spectra 
similar to our tadpoles. There were no systematic dif- 
ferences within 0.2 dex for 12 + log(0/H) > 7.7, which 
secures the reliability of the abundances found in most 
locations. As for the points with 12 + log(0/H) ~ 7.5, 
N2 overestimates the oxygen abundance, which again se- 
cures the low metallicity values we find. Since N2 pro- 
vides the O metallicity based on the flux of a N line, 
the metallicity may be biased in objects with unusual 
N/0. The effect of varying N/0 seems to be unimpor- 
tant too. Figure [11] shows the ratio as derived from 
[SII]A6717,6731/[NII]A6583 - Eq. (HH). (The [SIX] lines 
arc not available in two targets and therefore we can- 
not estimate their N content.) Even though the error 
bars are large, the observed N/0 looks rather constant 
along the galaxies, at approximately the plateau typi- 
cal of low metallicity galaxies (log(N/0) ~ —1.5; e.g., 
iPerez-Montero fc Continil |2009( ) . Obviously, if N/0 is 
constant then it cannot fake the observed metallicity 
drop. However, since the error bars of N/0 are large, we 
decided to run a chi-squared test (e.g.. lPress et aLlfigSGf) 
to further discard N/0 as the source of the measured O 
variations in KIS05639 and KIS06877. We compared the 
observed O with the fake variations to be expected if O 
is constant but N2 varies as N/0 in Fig. [TT] The test 
shows how the two variations arc inconsistent with 90% 
confidence. 

Figure 112b . shows the mass- metallicity relationship for 
our galaxies, plotting the oxygen abundance of the head 
versus the photometric mass of the galaxy. The most 
massive tadpoles tend to have the largest metallicities, 
although the metallicity values are displaced downward 
with respect to the mass-metallicity relationship found 
in the local universe. The dashed line shows the divide 
above which 97.5 % of t he loc al galaxies are found accord- 
ing to iTremonti et al.l (|2004l ). The tadpoles appear just 
below this line, so they are metal poor for their massefl 
The mass-metallicity relationship varies with redshift so 
that the same galaxy mass corresponds to lower oxygen 
metallicities at highe r redshifts (a 0.5 6 dex drop at red- 
shift 2.3, according to lErb et al.ll2006l ). Our tadpoles are 
on the mass-metallicity relationship observed at higher 
redshifts (see the dotted line in Fig. [T2k). 

We note that the metallicity of the heads also scales 
with their dynamical masses (Fig.[I2D). The relationship 
is even tighter than the relationship with galaxy photo- 
metric mass (Fig. 112b ). but we do not have enough points 
to judge whether or not the improvement is statistically 
significant. From a practical point of view, the reduction 
of scatter associated with the use of dynamical masses 
supports the reliability our estimate of this physical pa- 
rameter. 

8. NOTES ON INDIVIDUAL GALAXIES 

The general properties of the sample arc discussed in 
detail in the preceding sections. Here we focus on a few 
specifics of the individual galaxies. 

® As we already pointed out, the metallicity scale is not free 
from uncertainties. However, even taking them into account, the 
tadpoles are metal poor. The triple-dotted dashed line in Fig. 112b 
shows the solid line transformed to our metallicity scale using the 
prescription bv lKewlev fc EUisonl 1)20081 ). The observed metallici- 
ties arc well below this line too. 



KIS03193. This is the only tadpole that shows no rota- 
tion, and its line-width versus position curve has a cu- 
rious double-hump shape (Fig. [5l). Its light distribution 
is less lopsided than for the other tadpoles (Figs. [1] and 
|4|). The absence of rotation combined with the increase 
of linewidth toward the outskirts and the mild lopsided- 
ness may be consistent with a face-on disk, and in this 
sense KIS03193 differs from the rest of the sample. The 
observed linewidths are fairly large, comparable to those 
of other larger galaxies in the sample. Consequently, the 
dynamical mass of the head inferred from linewidths is 
significantly larger than its stellar (photometric) mass. 

KIS03867. In both mass and size, this is the smallest 
tadpole in the sample (see the photometric masses of the 
full galaxies in Table [2l and also the Ikpc scales given 
in Fig. |4]). Its linewidths do not differ so much from the 
linewidths of other tadpoles, which contrasts with the 
low photometric mass inferred in Paper I for the head of 
the tadpole. Consequently, the ratio of dynamical mass 
to photometric mass is particularly high in this case - the 
dynamical mass of the head turns out to be a hundred 
times larger than the stellar mass derived from photom- 
etry (Table [2]). The comparatively large oxygen abun- 
dance we derive for the head (12-1- log(0/H) ~ 8) is 
probably an overestimate due to observational errors. A 
few pixels away from the head the abundance drops by 
0.3 dex (Fig. [TUl), which would bring the galaxy down to 
a more natural location in the mass-metallicity relation- 
ship (Fig. [Hi). 

KIS05149. The asymmetry of its RC is remarkable: see 
Fig. [3] It extends much further out to one side, which 
would imply that half of the disk is missing from ob- 
servation. We cannot explain how this happens, unless 
the galaxy has a cigar-likc shape rather than an axi- 
symmetric structure. It is unlikely that large reddening 
obscures half the disk, since the SDSS image does not 
show a large color gradient across the galaxy. It may 
be a merger, but then it would have to produce an un- 
likely large scale velocity field resembling a RC. The oxy- 
gen abundance is fairly constant along the galaxy, which 
seems to be associated with being a massive object (it has 
the largest mass of the sample - Table [2]). The galaxy 
presents a second bright knot separated from the rotation 
center. 

KIS05639. It is one of the XMP galaxies in the sample, 
with the head being a young starburst (a few Myr old). 
It also presents a very irregular RC, particularly in the 
outer parts. 

KIS06669. This galaxy has a second bright knot, not 
far from the head (Fig. |4]). The two knots are of low 
metallicity, which contrasts with the metallicity of the 
rest of the galaxy (Fig. [T0| . The center of the RC stays 
in between the two emission peaks of the galaxy (and so, 
off-centered; see Fig. [3]). The linewidth curve shows a 
significant dip which does not coincide with any of the 
two emission peaks (Fig. [5]). 

KIS06877. This tadpole has the lowest metallicity and 
is the youngest object in the sample, and its bright head 
does not seem to contain dark matter. It also presents 
unusual high-excitation emission lines in the spectral re- 
gion of the WR bumps (Sect. [5]) . As wc explain in Sect. IH 
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Fig. 10. — Oxygen abundance variation across the galaxies. The vertical solid line represents the center of rotation, whereas the vertical 
dotted lines indicat e the location of maxima in the Hq flux profile. The thick horizontal solid line indicates the solar metallicity (from 
lAsplund et al.|[2009l ). Note the existence of abundance variations, with the minima coinciding with the largest Ha signals. Note also that 
KIS05639 and KIS06877 reach very low abundances, below one-tenth of the solar value, therefore, they are members of the set of rare XMP 
galaxies. The sorting of the galaxies and the abscissae are identical to those in Fig. \3\ 



we sought distortions in the RCs that could be associated 
with massive self-gravitating heads. Perhaps KIS06877 
represents the only example. Although we cannot as- 
sess the reliability of the distortion present in its RC 
at the head position (see Fig.|31), we carried out the aca- 
demic exercise of reproducing the observed velocity curve 
with two RCs combined, one for the galaxy plus one for 
the head. The center of rotation of the head has been 
forced to be given by the position of the head. The com- 
bined RC is shown as the dashed line in Fig. [T31 and it 
improves the fit to the data points (compare with the 
solid line). The head is modeled with a counter- rotating 
disk with a maximum velocity of 6kms~^, which ac- 
cording to Eq. ([5]) corresponds to a dynamical mass of 
log(A//AfQ) ~ 6.7. The mass thus derived is similar to 
both the photometric mass of the head and the dynami- 
cal mass inferred from lincwidth (see Table [2]). We also 
note the coincidence of this wiggle in the RC with an 
obvious decrease in line-width associated with the head 
(see Fig.©. 

KIS08466. The galaxy has a multi hump structure 
(Fig. [1] and [4]). The line-width curve presents a curious 
two-hump shape, which is not centered in photometric 



maxima (Fig. [5|). The oxygen metallicity is rather uni- 
form along the galaxy, as for the other massive galaxy in 
the sample, KIS05149. 

9. DISCUSSION AND CONCLUSIONS 

Galaxies with a bright peripheral clump on a fainter 
tail are called cometary or tadpole (Fig. [T]). The origin 
of the shape is unknown, but it may trace a transit phase 
in the assembly of many disk galaxies (Sect. [T]). Low 
mass local tadpoles were identified and studied photo- 
metrically in Paper I. Here we follow up the study, and 
analyze the chemical and dynamical properties of seven 
such targets inferred from long-slit spectra around Ha. 

Five out the seven observed tadpoles show evidence for 
rotation 70%), and a sixth target hints at it. Often 
the center of rotation is spatially offset with respect to 
the tadpole head (three out of five cases) . The RCs of the 
smaller targets are not smooth but present fluctuations, 
suggesting a complex dynamics (e.g., a counter-rotating 
head - see Fig. [13]). The size and velocity dispersion of 
the heads are typical of giant HII regions and follow the 
scaling relationship known to exist between these two 
quantities (Fig. [6]). We find changes of velocity disper- 
sion along the galaxies, but they are not correlated with 
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Fig. 11. — Variation of the ratio N to O across the galaxies. The horizontal sohd hne represents the plateau at some log(N/0) ~ —1.5 
found in most low mctallicity galaxies, for which N is produced as a primary element. The sorting of the galaxies and the abscissae are 
identical to those in Fig.|3l The ratio could not be estimated in the two galaxies of largest redshift, where the required emission hues lie 
outside the observed spectral range. 



intensity variations. If anything, there is a tendency for 
the tadpole heads to coincide with local minima of veloc- 
ity dispersion. The head is defined as the position on the 
galaxy with the largest surface brightness. Observation- 
ally, we find it to coincide with the region of largest Ha 
flux, and so, of largest SFR in the galaxy fe.g.. lKennicuttl 
Il998f ). Moreover, we also find the continuum fiux to ex- 
tend further out as compared to the flux in Ha, which 
is concentrated around the head. Thus the bright heads 
seem to be large starbursts with their random motions 
reduced with respect to the rest of the galaxy. 

Using the observed RCs and velocity dispersions, we 
estimate the dynamical masses of the galaxies and their 
heads. The dynamical masses of the full galaxies are 
between three and ten times larger than the stellar mass 
inferred from photometry. The dynamical masses of the 
heads also exceed their stellar masses, but to a lesser 
extent than the full galaxies. Actually, the photometric 
mass and the dynamical mass of three heads agree within 
error bars. In two other cases, however, the dynamical 
mass of the head exceeds the photometric mass of the 
full galaxy. 

The oxygen metallicity estimated from [NII]6583/Ha 
often shows significant spatial gradients across the galax- 



ies (~0.5 dex), being lowest at the head and increasing 
in the rest of the galaxy, tail included. So far as we are 
aware of, this is the first time that a metallicity growing 
away from HII regions has been reported in local galaxies. 
The sense of the resulting metallicity gradient is at vari- 
ance with the observation of local disk galaxies, where the 
gas-phase metallicity increases toward th e galaxy centers 
(IVilchez et allllQSSt iGarnett et al.l[l997f ) or is just con- 
stant ( Moran et al.ll2012h . However, the type of variation 
we measure, with a minimum metallicity at the most in- 
tense star-forming regi on, has been observed in galaxies 
at redshift around 3 bv lCresci et al.l (|20100 where it is in- 
terpreted as evidence for infall of pristine gas triggering 
star formation. Once systematic errors are disregarded 
(as we did in Sect. [7]), it is difficult to avoid such interpre- 
tation, also in the case of our tadpole galaxies. We con- 
sidered and then discarded the following two possibilities: 
(1) assume a regular metallicity distribution decreasing 
outward. If the head was formed from gas in the galaxy 
outskirts, but has spi raled in toward the cen ter by dy- 
namical friction (e.g.. lElmegreen et al.l[2012at ). then the 
head would naturally present a metallicity lower than its 
immediate surroundings and similar to the galaxy out- 
skirts. However, this prediction of metallicity gradients 
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Fig. 12. — (a) Mass metallicity-metallicity relationship for the 
galaxies. We consider the photometric mass of the galaxy and 
the oxygen abundance at the galaxy head. Th e mean relationship 
found in the local universe by riVemonti et al.l II2004I ') is shown as 
the solid and the dashed lines — the thick solid line represents the 
mean of the distribution, whereas the dashed line indicates the line 
above which 97.5 % of the local galaxies arc found. Our targets are 
outliers of the local relationship; too metal poor for their masses. 
They are more in the line of t he starb urst galaxies found a high 
redshift (the dotted line, from lErb et a l. 2006). The triple-dotted 
dashed line corresponds to the solid line transformed to our m etal- 
licity scale using the prescription bv lKewlev &c EllisonI I I2008I ). (b) 
Same as (a) using the dynamical mass of the heads. The solid 
and dashed lines arc also the same, and have been included for 
reference. 
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Fig. 13. — Velocity curve of KIS06877 fitted by a single RC (solid 
line) and including a counter-rotating disk at the position of the 
head (dashed liric). The data points and the solid line are identical 
to those in Fig. O 



induced by internal migration is inconsistent with our 
observation, where the metallicity at the head is not just 
lower than the surroundings but the lowest (see Fig. [TO)) . 
(2) Metal-rich supernova (SN) driven winds remove met- 
als from shallow gravitation al potential wells, produc- 
ing metal-pqor ga laxies (e.g.. iMac Low fc Ferraral[l999l : 
iRecchi et al.l 1200^ . This mechanism explains the over- 
all low gas metallicity of some galaxies with significant 
old stellar populations, but it does not account for the 
presence of a region like the tadpole head, with metal- 
licity lower than the rest of the galactic gas. The metals 
ejected by the winds are those created by the stars that 
explode as SNe. They are not the metals of the ambient 
gas. Then the loss of these ejecta reduces the efficiency 
of the galaxy to retain metals, but it does not reduce 
the metallicity of a particular region of the galaxy. In 
short, the interstellar medium of the tadpole galaxies is 
not well mixed but shows significant metallicity gradi- 
ents. Since the mixing time scale is expec ted to be rela- 
tively short (shorter than a few My r; e.g.. lTenorio-TagT3 
Il99a Ide Avillez fc Mac Low! 120021) . the recent infall of 
metal-poor gas seems to be the only viable alternative 
to explain the metallicity drop observed at the tadpole 
heads. 

The observationally motivated interpretation of exter- 
nal gas infall fits in well the cold-flow gas accretion sce- 
nario arisiiig_Jrom cosmolo gical numeri cal simulations 
(e.g.. lKeres eTa l. 2005; Dck cl et al.ll2009[ ). It predicts lo- 
calized accretion in clumpy streams of pristine gas ready 
to create stars. The streams may directly form giant 
clumps that we detect as tadpole heads or, alternatively, 
feed the disks with turbulent gas that eventually frag- 
ments into giant clumps by gravitational instability. In 
both cases the massive clumps are prone to migrate 
toward the galax y centers and become progenitors of 
centr al spheroids (jElmegreen et al.ll2008l : iCeverino et al.l 
l2010f ). Other details of the observed tadpole properties 
are also consistent with the cold-flow accretion scenario. 
The process is expected to be ubiquitous at high red- 
shift. Then the flows fade away gradually in a process 
that has not bein g completed in small galaxies yet (e.g., 
iKeres et al.ll2009| ) . This prediction is consistent with the 
absence of low metallicity HII regions in large local spi- 
rals, as well as in our most massive tadpoles (kiso5149 
and KIS08466; see Fig. [10]), for which the infaU would 
be already over. The metallicity drop is observed only in 
the low mass objects, reflecting the downsizing process 
in galaxy formation. 

The geometrical displacement of some of the tadpole 
heads with respect to the centers of rotation also favors 
the cold-flow scenario. The expected streams of cold 
gas never end up at the galaxy center. The clumps are 
formed in the disk, and require time to be transported 
inward. 

Extremely metal poor (XMP) galaxies are rare. Tad- 
pole galaxies are also rare. The fact that we observe two 
XMP galaxies in a sample of seven tadpoles cannot be 
casual (Sect. [T0| . It is known that a signiflcant fraction 
of XMP galaxies turns out to be tadpole or co metary 
(jPapaderos et al.l 120081 iMorales-Luis et al.ll2011| ). Here 
we flnd the reverse to be true as well, i.e., tadpoles have 
a significant chance of being XMP. The coincidence of 
these two seemingly disconnected properties is best un- 
derstood if the objects are primitive, with the cometary 



16 



shape and the low metalUcity reflcxing dynamical and 
chemical youth, respectively. 

All these results combined are consistent with the local 
tadpole galaxies being turbulent disks in early stages of 
assembling. Their star formation seems to be sustained 
by accretion of external metal-poor gas. 
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